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Abstract. We discuss the crossover between the small and large field cutoff 
(denoted x ma x) limits of the perturbative coefficients for a simple integral 
and the anharmonic oscillator. We show that in the limit where the order 
k of the perturbative coefficient Ukixmax) becomes large and for x max in the 
crossover region, a^(xmax) oc J ma ' c e — A ( x ~ *o(*0) d x . The constant A and 
the function xo(k) are determined empirically and compared with exact (for the 
integral) and approximate (for the anharmonic oscillator) calculations. We discuss 
how this approach could be relevant for the question of interpolation between 
renormalization group fixed points. 
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1. Introduction 

The existence of different behaviors at different scales is a central difficulty in quantum 
held theory. This is clearly the case for QCD [I| where the short distance behavior can 
be described using perturbation theory while confinement and other non-perturbative 
phenomena appear at large distance. A similar situation is encountered in scalar field 
theory and spin models where, for instance, descriptions of the renormalization group 
(RG) flows near the high-temperature fixed point in terms of perturbative expansions 
about the Gaussian fixed point are in general problematic. 

In general, the fact that aspects of the RG flows cannot be described using 
perturbation theory is not surprising given that perturbative series usually have a zero 
radius of convergence [21 ■ The convergence of perturbative series change drastically 
[13 E] if the large field configurations are removed from the path integral. This can be 
implemented, for instance, by removing configurations such that the norm of local 
fields exceeds some cutoff value. The field cutoff modifies the original problem, 
however, it has been found that for nontrivial <f> 4 problems 0] the modified series 
apparently converge toward values which are exponentially close to the exact ones. In 
addition, at fixed order, it is possible in simple examples to choose the field cutoff in 
order to optimize the accuracy EJ- 

The calculation of the perturbative series for the modified theories where a large 
field cutoff is introduced is non-trivial. In the case of the anharmonic oscillator, 
approximate analytical methods [7] have been developed in the limits of large and small 
field cutoffs. Except at very low order, these methods are not accurate in the crossover 
region between these two regimes and it is imperative to develop new methods that are 
accurate in the crossover region. Fortunately, remarkable regularities were observed 
empirically: with appropriate rescalings and translations, the functions expressing the 
perturbative coefficients in terms of the field cutoff approximately collapse along a 
single curve. This suggests that despite the complexity of the perturbative coefficients, 
a simple description of the crossover seems possible. 

In this article, we refine the description of the approximate collapse of the 
crossover functions discussed above. An important feature is that the collapse 
improves as the order increases and a universal (order independent) function is reached 
asymptotically. In order to model properly this approach, we start with a simpler 
example, namely the perturbative series associated with a one variable integral. This 
is done in Sec. El where a very good agreement between an empirical polynomial 
parametrization and a saddle point approximation is found. Generalizations of these 
results are discussed in Sec. |3J Our main focus in this article is to discuss the case 
of the anharmonic oscillator "with a field cutoff". This model is introduced in Sec. 01 
and recent results regarding the behavior of the perturbative coefficients 7J are briefly 
reviewed . Approximate Gaussian parametrizations of the crossover are introduced in 
Sec. Higher orders corrections are discussed in Sec. H3 where a simple description 
of the crossover for large order emerges. The possible relevance of this approach for 
the description of RG flows is discussed in the conclusions. 

2. The simple integral 

In this section, we discuss the calculation of the simple integral 




(1) 
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In this example, introducing a field cutoff amounts to remove the integration tails 
where \x\ > x max . Their contributions to the integral is less than y/2Tre~ XXmam and if 
Xmax is large enough, we can approximate Z(X) with 

Z(X,x max ) = f ^ dxe-^ 2 - Xxi . (2) 

If x max is kept at a fixed finite value, the Taylor expansion of the exponential 
converges absolutely and uniformly over the domain of integration and it is legitimate 
to interchange the sum and the integral. The resulting series converges over the entire 
A complex plane and reads: 

oo 

Z(X,x max ) = y^Qfc(x max )A fc , (3) 

fc=0 

with 



(-1)' 



a k (x max ) = ^— p- I dxe * x x Ak . (4) 



On the other hand, in the limit where x max — > oo, we have 

a fc M = ^r(2fc + i)(2) 2fe +^ , (5) 

and the coefficients grow factorially with the order which implies that the series 
diverges for any non-zero value of A. In Fig. ^ we display the ratios 

Rk {Xmax ) = 7 \ (6) 

a k (oo) 

for k = 1 to k — 10. As we can see, the curves corresponding to the different orders 
have a similar shape and move to the right as the order increases. Each curve can 
be characterized by three regimes: the low field cutoff regime, where the coefficient 
becomes very small, the large field cutoff regime where the coefficient reaches its 
asymptotic value and the crossover region. We now discuss each case separately. 

2.1. The a; max — > oo limit 

When x max — * oo, we can use integration by part recursively for the missing tails of 
integration in Eq. QJ: 



1 - Rk{x max ) ~ e ^ 



1 .,,2 



r/1 2 \2k-± (12 \2k-i 
\ 2 max/ . V2 max/ 



(7) 



r(2fc + i) r(2fc-±) 

At each integration by part, the argument of the T function at the denominator 
decreases by one. Since the arguments are half integers, the process does not terminate. 
Using l/r(x — 1) = (x— 1)/T(x), we end up with a factorial growth for the coefficients 
of the power series. 

2.2. The a; ma x limit 

In the opposite limit where x max — > 0, we expand the exponential in Eq. Q and we 
obtain the converging series in x^ax '■ 

H k (X max) - r(2fc+ l ) Z^ n(4fc + 2Z+1) • W 
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Figure 1. The ratio Rh{xmax) as defined in Eq. JSJ. As k goes from 1 to 10, 
the curves move to the right. 




Figure 2. The ratio Rk(xmax) in the large x max and small x max approximations. 
The dot stands for the accurate numerical result. The dash lines correspond 
to the large x max approximation. The solid lines correspond to small a; max 
approximation. The labels 1, 2, 3 correspond to the number of terms in Eq. 
171 . The six lines on the bottom left correspond to polynomial expansions when 
£max — > from Eq. JgJ with truncations at order 0, 5, 10, 15, 20, 25. 



It is clear that the factorial growth of the denominator implies that this series converges 
over the entire x max complex plane. 

In Fig. [21 we show the curves corresponding to different orders in the two 
expansions discussed above, for k = 5. The small x max region needs to be resolved 
logarithmically. If we restrict the expansion in Eq. to I < L. for some given L, 
we notice that the truncated series becomes a poor approximation when x max reaches 
some critical value (almost independent of the order k in A). This is illustrated in Fig. 
121 for L = 10, where the critical value is near x max ~ 2.7. One can roughly estimate 
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this critical value by requiring that the L-th term is the same order as the L — 1-th 
term. This yield the critical value x max ~ y2L. A more accurate way to proceed is 
to determine the value of x ma x for which the L-th order approximation becomes poor. 
The critical values of x max can be fitted reasonably well with 0.63 + 0.64VL. We can 
thus conclude that despite the fact that the series is convergent, many terms need to 
be calculated when the order becomes large. The number of terms that we need to 
calculate in order to get an accurate answer for a given x max grows like x\ 



max ' 



2.3. Approximate universality 

Fig. and the above discussion show that it is not easy to use the two expansions 
in the crossover region, which corresponds to the region where the integrand of Eq. 
JIJ peaks. Using a saddle point approximation about the maximum of the integrand 
x = v4fc, we obtain 



J"*""- dxe~ 2k+4kln ^ 



Rk{x m ax)^ joo dxe - 2k +ikHVik) e -^-Vik) 2 + - ' ^ 

If we only retain the quadratic term in the above expansion of the phase, we obtain 
what we call hereafter the Gaussian approximation. This approximation fails to 
reproduce the small x max behavior of Eq. |JSJ. It is clear that far from the peak 
of the integrand, the Gaussian approximation is not good but it is in regions where 
the original integrand is very small. One can simplify this expression by adding the 
left tails of integration. This makes sense for k not too small, where it generates errors 
of the order (4\/fc)~ 1 e _4fc . For sufficiently large k, the following approximation holds 
for x m ax not to far from \/4fc 

Rk{x ma x) * 7T- 1 / 2 f ^ dxe-^^ . (10) 
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In Eq. i|10|) . the only dependence on k is in the argument of the exponential. Its 
net effect is a translation in x max and the shape of the function is universal (k- 
independent). In other words, 



) =s U int { V4fc) (11) 

with 

U mt (x) = J -°° " . (12) 

If the assumption of universality expressed in Eq. (|11J) is approximately correct, the 
data should collapse into the universal curve when the argument of the fc-th order 
curve is shifted by V4fc. Fig. 0] confirms that this is approximately the case. A more 
detailed look at Fig. 0] indicates that as the order increases, the shifted curve gets 
closer and closer to Umt- This limit is studied in the next section. 

2.4- Order- dependent corrections 

The basic approximation used to obtain Ui n t in the previous section is the saddle point 
approximation of the integral. We expanded the phase of the integrand 

^(w) = hxRfc(aWr;) (13) 

about its maximum up to second order. It is possible to improve this approximation 
by expanding to higher order 



-x 2 +4k\nx = 2k(21n(2Vk)-l)-(x-VIkf-4:k^2(l/r 

71— 3 



(V4fc-x)/V4fc 



.(14) 



The series converges for < x < 2 v 4k. The even terms are negative and if we truncate 
at even orders, the integral of the exponential of the truncated sum converges when 
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Figure 5. Values of the third and fourth coefficients of the Taylor expansion of 
the phase about its maximum obtained from polynomial fits compared with the 
predictions of Eq. 1141 



the bounds of integration are sent to infinity. As the order increases, the contribution 
of the two regions outside the radius of convergence are effectively cutoff, because the 
phase $fc becomes large and negative in these regions. For x < 0, this is desirable 
because there is no contribution from this region in the original integral. On the other 
hand, contributions from x > 2y4k exist in the original integral. However, as we move 
from the maximum of the phase at \/4Jc to the edge of the radius of convergence at 
2\/4fc, the phase of the integrand drops by (— 6 + 41n2)fc ~ — 3.2k and the contribution 
of the right tail in the original integral is small, even for relatively small values of k. 

In order to check the validity of this approach, we have calculated the phase 
by using a discretized approximation of the derivative of Rk- We then fitted a set 
of points near the maximum with polynomials. The results were quite robust under 
changes in the number of points used, the degree of the polynomial and the spacing 
used to calculate the derivative. The second coefficient stabilized rapidly between 
-0.99 and -1 in agreement with Eq. (|14|l . The coefficients of order 3 and 4 are shown 
in Fig. [S] where the good agreement with Eq. I|14l) is quite obvious. 

The inclusion of the corrections of order 3 and 4 in the phase significantly reduces 
the small differences between Ui n t an d the data in Fig. 21 and, for instance for k = 5, 
it is impossible to distinguish between the corrected integral and the data with the 
naked eye. 
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It is now clear how and why the agreement between the data and the simple guess 
U int of Eq. ijT2|) gets better as the order increases. The size of the crossover region 
is of order 1 and controlled by the quadratic term in the phase. In this region, for 
\x — \/4fc| ~ 1, 3?fc — const. — (x — V4fc) 2 because corrections of order n > 2 in Eq. 
<|14fl are suppressed by a factor fc 1- "/ 2 . On the other hand the argument does not 
hold far away from the crossover, for \x — \^4k\ ~ Vk, where a large number of terms 
is necessary to approximate logarithms in the phase that reproduce the correct power 
behavior of the integral. 



3. Remarks about the saddle point approximation 



In this section, we generalize some aspects of the saddle point approximation used in 
the previous section. One important feature of the expansion of the phase given in 
Eq. (|14|l is that the coefficient of the quadratic term is fc-independent and twice the 
value of the coefficient of the quadratic term in the original integral. This has a very 
simple explanation. In general, e~ Ax x B is maximal at x* = \J ' B/2A and we have the 
expansion 

e~ Ax 'x B cx e- 2A< - x ~ x *) 2+ - . (15) 

In doing this calculation, one realizes that in the second derivative of the log 
estimated at x* , the factors B cancel exactly. We want to emphasize that this "non- 
renormalization" is due to the fact that the exponential is multiplied by a single power. 
If we consider instead e~ Ax (x B + ex B+s ), then the exact cancellation does not hold. 
This can be checked at first order in e, where x* is changed according to 



y/B/(2Aj - VB/(2Aj 



l + (e5/(2B))(B/(2A)) s ^ 



(16) 



and the coefficient of the quadratic term in the expansion of the phase is changed 
according to 



2A 



-2A 



1 - (eS 2 /{2B))(B/(2A)) s ^ 2 



(17) 



and becomes independent (except if 8 = 2 or 0). 



4. The anharmonic oscillator with a field cutoff 



Our main objective in this article is to understand the crossover for the anharmonic 
oscillator. The model has been studied in detail in Ref. |7]. We introduce some 
notations and briefly review some basic results. For convenience, we use quantum 
mechanical notations instead of field theoretical ones. The "field" will be denoted x 
instead of <fr and the field cutoff will be denoted x max . In quantum mechanics language, 
it means that the potential becomes infinite at ±x max - We use units such that h, lo 
and the "mechanical mass" m are 1. If the dependence on these 3 quantities were 
restored, dimensionless quantities would be expressed in terms of ^Lum/hx max and 
hX/m 2 uj 3 . The Hamiltonian reads 

H=£+V{x) (18) 

with 
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We will consider the perturbative expansion of the ground state: 

oo 

E {x max X) = $^ fc) (a: m(M )A fe . (20) 

k=0 

As in section [2J we use the notation Rk for the fc-th coefficient in units of its usual 
value: 

Rkipmax) = E^(x max )/E^ k \oo) . (21) 

We have shown |7j that in the limit of small x max , 

Rk(Xmax) °C X max , (22) 

while in the limit of large x max , 

1 - i? fc (x max ) oc x^e-^ . (23) 

5. Gaussian approximations of the crossover 

In this section, we discuss Gaussian approximations of the crossover in a way similar 
to what was done in Sec. 12.31 for the integral. From Eq. (|23|) . we expect that in the 
large x max limit, we have at leading order 

Rhi. x max) OC ^-max e max • (24) 

If we now extrapolate backward this behavior at intermediate values of x max (this is 
of course an approximation), R' k (x max ) has a maximum at 



x*(k) = V2k + 1 . (25) 

Expanding the logarithm up to order 2 about this point, we obtain the Gaussian 
approximation: 

R' k {x max ) cx e -2(x max -Mfc)) 2 . (26) 

With this approximation, we have R' k '(x*(k)) = 0. This approximate prediction 
for the value of x for which the second derivative vanishes can be compared with the 
accurate numerical value, denoted xo(k), for which it actually vanishes. The numerical 
values can be found in Table 2 of Ref. and are shown in Fig. EJ together with 
approximate analytical formulas. Typically, xo(k) — x*(k) with a 10 percent accuracy 
for k < 10. A more accurate formula can be found by fitting xo(k) 2 with a linear 
function. Using the data for k = 2 to 10, we obtain Xo(k) ~ Vl-63fc + 3.26. Using the 
data for k — 10 to 20, we obtain a slightly different fit that we expect to be closer to 
the asymptotic behavior: 



x (fc) ~ Vl.61fc + 3.48 . (27) 

In Ref. |7], we found that if Rk{x max ) is translated by xo(k), the data for 
k = 2, . . . 10 approximately collapses. In other words, 

R k (x + x (k))~U anh (x) . (28) 

An approximate form for Uanh{x max ) is suggested by Eq. (|26|l . Integrating and 
normalizing by requiring lim x ^ +co U anh {x) - 




U anhtl (x) ~ / dye-* . (29) 
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Figure 6. The numerical values of xo(k) (circles) compared to yJ2k + 1 (dash 
line) and V1-63A; + 3.26 (solid line). 



We have used the subscript anh, 1 to indicate that it is a first approximation and that 
a better form will be found later. This proposal is compared with the empirical data 
for k — 7, . . . , 10 in Fig. [7| We only displayed a few values of k so that it is possible 
to distinguish the various orders on the graph. One can see that the collapse of the 
various curves is reasonably good. The function U an h,i fits the collapsed data quite 
well on the upper right, but appears to be slightly above the data on the lower left. 
However, as the order increases (in the Figure from 7 to 10), the difference diminishes. 

The reason why y/2k + 1 is not a very good approximation of Xo(k) is that the 
large x max behavior given by Eq. (|24J) is not accurate for intermediate values of x max . 
Studying numerically the behavior of 

Rk{x m ax)e x ™°* , (30) 

near xo(k), ones obtains a very good linear behavior on a log-log plot. If Eq. (|24|l were 
correct, the slope would be 4fc + 2. However, collecting the slopes for various values 
of k, one concludes that 4fc + 2 should be replaced by 3.2k + 7.8 in Eq. (|2"4"jl . This 
is consistent with the empirical form of xo(k) of Eq. I|27|l which is approximately the 
square root of one half of this fit (see Sec. |3]with A—l.) When studying the difference 
between Eq. (|30() and its approximation by a single power, one notices small quadratic 
corrections. We will not report all the details of this investigation, but only mention 
that these quadratic corrections can be approximately removed by replacing Eq. (|30[) 
by R' k (x max )e^ kXma:c and adjusting the value of (3 k to the value for which the quadratic 
corrections change sign. 

We conclude that a more realistic version of Eq. I|24|) is 

Rhi. x max) OC X max ^ e k max , (31) 

The results of Sec. suggest that e~ 2y should be replaced by e~ 2l3kV in Eq. 
(|29|l . This brings an explicit dependence on k, and universality can only be reached 
asymptotically if /3k tends to a limit when k becomes large. 
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Figure 7. Ri c (x^.xo(k)) for k = 7, ...10 and the function C/ an j 1 i(a;) of Eq. 1291 , 



Empirically, the values of (3k are slightly larger than 1 for k < 10 and slightly 
smaller than 1 for k > 10. This modification allows to reduce the small discrepancies 
of the Gaussian approximation used before. Note that in Fig. [7J the closeness of k = 9 
and 10 to U an h,i reflects the fact that (3g and /3io are very close to 1. (3k is closely 
related to /2, a quantity that will be studied in the next section. Anticipating the 
results presented there, we expect that as k increases, (3k stabilizes to a value close to 
0.75. 

The limitations of the Gaussian approximation, or of any other expression of a 
universal function U as the integral of a even function, can be seen from the duality 
relation which exchanges the small and large field cutoff regions 

U(-x) = 1 - U(x) , (32) 

which cannot be exact since the approach of when x max — > is power like while in 
the large x max limit, the approach of 1 is exponentially small. In at least one of these 
two limits, U should fail to provide the correct behavior. 

6. Numerical study of the higher order corrections 

The complexity of the exact form of the terms of the perturbative series increases 
rapidly with the order (see for instance |S] for exact expression up to order 4). The 
number of graphs grows factorially 9 . It seems unlikely, that at this point, one 
could guess approximate equations such as Eq. (|31J) . or an improved version of it, 
from analytical expressions at large order. For this reason, we will describe \n(R' k ) 
near xo(k) using polynomial approximations. This method is numerically robust and 
has been tested successfully in the case of the integral where the expansion can be 
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calculated analytically. Our main working assumption is that in the vicinity of x a (k): 

Q 

ln(4(x)^4 9) (,-x„(fe))' (33) 

q=0 

Note that the presence of odd powers in the expansion breaks the duality expressed 
in Eq. i|32|) . The fitted value for the largest order (q = Q) coefficient is typically not 
very stable and as we want to focus on the terms with q = 2, 3 and 4, we have used 
Q = 5 and 6. We then checked that the terms A k for q = 2, 3 and 4 were quite 
stable under small changes as in Sec. El They are displayed in Fig. |SJ 

As the A k 9 ^ seem to reach asymptotic values, we have fitted the larger fc values 
with the parametric form 

A^=a^+b^k-° M (34) 

Other parametric forms such as an exponential decay have been tried too but have 
been shown to have a significantly larger chi-square. We first estimated by making 
a log-log plot of — A k q ^\ and then obtained and with standard linear 

fit methods. We used different sets of A k with k running from ko to 20. The results 
showed a regular dependence on kg: as fco increased, the three fitted parameters 
showed small, monotonous changes and we then attempted an extrapolation (to large 
fco) using the same method as above. This second extrapolation is of course much 
more delicate and used only to estimate the errors. The stability of the procedure was 
checked by using five different sets of A k 9 ' obtained with slightly different procedures: 
changes in the range of x max used in the fits, changes in the Ax max used to calculate 
the derivative of Rk, changes in Q = 5 or 6. 

For q=2, the fits at fixed fc were found to be quite consistent among the five sets 
and the second extrapolation was in some cases possible. Our final estimate is 

4 2) - -1.5(1) + 2.0(l)£r°- 6(1) , (35) 
with the parenthesis indicating the estimated error on the last written digit. For 
sufficiently large fc, we found that A k ~ —2(3k as defined in the previous section. It is 
clear form Fig. [SJthat for fc ~ 10, A k ~ —2 in agreement with the observation (3 W ~ 1. 

(2) 

To give an idea about the accuracy of the agreement, for fc = 15, A 15 ~ —1.866 while 
0i 5 = 0.938 = 1.876/2. 

For q = 3, more variations among the five data sets were observed. A tentative 
value for a*- 3 -* is 0.02(1). However this is more than 4 time less than A 1 ^ ~ 0.089(1). 
More stable results with consistent extrapolations were obtained by setting a^ = 0. 
Our final estimate is 

4 3) s 1.3(3)fc~ a9(1) . (36) 

We repeat that it is difficult to rule out a small positive value asymptotically. The 
decay rate seems faster than in the single power case studied in Sec. where 
A k 3 ^ oc fc -1 / 2 . In our fits the average value of the power was -0.94 for fco = 17 
and the extrapolated value the farther away from this average was -0.86. 

For q = 4, stronger variations among the five sets were observed. In all cases a*- 4 -* 
was much smaller than A^q — —0.017(1) and with no consistent sign. More stable 
results with reasonably consistent extrapolations were obtained by setting = 0. 
Our final estimate is 

4 4) ~ -1.5(6)fc- 1 - 3 ^ (37) 
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Figure 8. for q = 2, 3 and 4 and k = 1, ... 20. 



We have also considered the coefficients for larger values of q and found that they 
also appear to reach zero asymptotically. Their values are in general quite small. For 
instance, A$ ~ 0.004 and A$ ~ -0.0007. 

We conclude that for k sufficiently large and x max = xo(k) + x, with x of order 
1, we have in good approximation that 

R k (x {k) + x) ~ U anh (x) ~ V2l7^ ^ dye- 2 ^ 2 , (38) 
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with A ~ 0.75. The words of caution regarding this "limit" are the same as in the 
case of the integral: if |x| is too large, the non-Gaussian corrections need to be taken 
into account. 

7. Conclusions 

In summary, we have considered the simplest field theoretical models with a field 
cutoff, in zero and one dimension. We found in both cases, that the crossover behavior 
of the perturbative coefficient can, for orders large enough, be described in good 
approximation by rescaling and translating a universal (order-independent) function. 
Furthermore, this function can be expressed as a Gaussian integral. This statement is 
exact for the simple integral (in the sense of Eq. (|14[>). For the anharmonic oscillator, 
small non-Gaussian corrections persisting at large order cannot be ruled out. 

We expect similar features in higher dimensional scalar field theory with a UV 
regulator. Generically, a modified perturbative calculation involves a few coefficients 
with values close to their usual ones, a few coefficients in the crossover region, and 
the rest of the coefficients taking values much smaller than in regular perturbation 
theory. These three regimes are reminiscent of the three regimes encountered when 
calculating renormalization group flows between two fixed points. One should notice 
the similarities between the graphs presented here and those of Refs. ^3 where 
this type of crossover behavior is studied. In the same way, the approximate duality 
of Eq. il-il'li reminds of the duality between fixed points of a simplified renormalization 
group equation |12j or the exchange between the perturbative and non-perturbative 
sectors of quasi-exactly solvable periodic potentials ^3] ■ 

We plan to make these analogies more specific (or disprove these ideas) , by making 
explicit calculations with Dyson's hierarchical model. Namely, we plan to compare 
the relative weights of the three types of perturbative terms discussed above with the 
relative increases in the three regions of the renormalization group flows discussed 
above for the zero-momentum n-point functions. Our hope is that we will be able to 
trade the very complicated problem of the interpolation between RG fixed points to 
a more tractable one. 
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